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FOREWORD 


Holography appears to have the potential of providing a unique measurement 
procedure for use in cloud microphysical research because of its ability to capture and 
reproduce a three-dimensional image of an event without interfering with the event. The 
motion, velocity, and growth cycles of particles can be recorded in a manner not 
obtainable with standard photographic techniques used in most cloud physics research 
programs. The study reported herein resulted from meetings and discussions with the 
author and Mr. William W. Vaughan, Chief, Aerospace Environment Division, 
Aero-Astrodynamics Laboratory, Marshall Space Flight Center, and others engaged in 
feasibility studies for a zero-gravity (zero-g) atmospheric cloud physics laboratory as a 
payload for the Space Shuttle/Spacelab. The study documented by this report was made 
to investigate the feasibility of applying holographic techniques to the real-time 
determination of atmospheric cloud particles parameters in a zero-g environment. 

The Zero-g Atmospheric Cloud Physics Laboratory is being developed under the 
technical direction of the Marshall Space Flight Center on behalf of the NASA Office of 
Applications and the Office of Manned Space Flight. The primary objective of this 
program is to provide a multi-experiment laboratory which members of the cloud physics 
scientific community may use in their research of cloud microphysical processes. 
Emphasized will be processes for which a near zero-g environment will enhance their 
capability to improve our current knowledge for application to cloud physics and weather 
modification problems. 


William W. Vaughan 
Project Manager 
Zero-g Atmospheric Cloud 
Physics Laboratory Project 
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OPTICAL HOLOGRAPHY APPLICATIONS FOR THE ZERO-g 
ATMOSPHERIC CLOUD PHYSICS LABORATORY 

INTRODUCTION 


Statement of the General Problem 

The science of meteorology has advanced rapidly toward the development of an 
understanding of large-scale atmospheric processes. Large-scale flow patterns in the 
atmosphere can be predicted with reasonable success because of the progress in defining 
and refining the equations of motion and conservation of mass and momentum. These 
large-scale flow patterns enable one to define and predict patterns of cloudiness. This 
cloudiness itself, the processes in the clouds, and the redistribution and release of energy 
in the clouds are what the average citizen thinks of as weather. 

In nature there are many cubic kilometers of clouds. Most cloud elements are 
inside the cloud and are not affected by the edges of the cloud. A terrestrial 
experimenter has the problem of trying to duplicate a system in which there are no edge 
effects and in which a cloud can survive long enough to have a number of things happen 
to it. The smallest of these cloud particles fall at several centimeters per second, and as 
they grow they reach fall velocities above 1 m/sec. This means that in a terrestrial cloud 
chamber that is 1 m in the vertical dimension, there is less than 1 sec of observational 
time for the larger of these particles. Attempts to overcome this limitation include 
capturing individual particles and placing them on wax paper. Teflon, copper, or stainless 
steel surfaces; hanging them on a thermocouple; suspending them on a spider web or a 
thread; or placing them between two immiscible liquids. These approaches have not been 
very successful because the suspension medium generally causes effects greater than the 
forces or actions that are being measured. 

Many significant accomplishments have been made in terrestrial laboratories. 
However, a comparison between the results achieved in the laboratory and what is 
observed in nature often gives no correlation whatsoever. In studying the whole regime of 
meteorology from the large-scale motions which produce cloud systems down to what is 
called weather (i.e., rain, snow, lightning), there are gaps in understanding the 
microphysical processes that occur between the inception of the cloud system formation 
and the events eventually occurring at the ground. Many of these data lie in the area 
called cloud physics. Meteorologists have long understood the requirement for a much 
better understanding of cloud microphysics. For the last 30 yr concentrated efforts have 
been made to understand some of the following atmospheric processes: Why does one 
cloud develop a spectrum of broad droplet sizes while others develop narrow spectra of 
sizes? Why does one cloud precipitate and others do not? Why does one cloud develop 
rapid electrical charging, charge separation, lightning and thunder, while other clouds 



similar in outward appearance do not? Why does one thunderstorm produce hail and 
another not? Why does a field of thousands of clouds produce hundreds of moderate 
thunderstorms, but only one develops into a tornado? Microphysicists, physical chemists, 
and applied physicists, as well as meteorologists, are still heavily involved in the study of 
the phenomena associated with these atmospheric processes. 

The range of cloud physics laboratory research extends from the millimeter rain 
drops and ice crystals to submicrometer condensation nuclei. Nature requires at least a 
million 10-pm-diameter cloud drops to combine in order to produce a 1 -mm-diameter 
precipitation drop. These droplets possess certain surface, electrical, and aerodynamic 
properties which establish whether these million small drops can combine to form the 
one big drop. Individual ice crystals must also be studied to determine how they grow 
and whether these ice crystals will somehow splinter and multiply to form more ice 
crystals. Other problems include electrical charging of ice crystals during the growth and 
collision processes and the effect of this charging on thunderstorm electrification. Many 
microphysical processes have been studied in some detail in the laboratory but under 
conditions which are not very representative of those in a cloud. In any cloud, the 
average element or group of elements within that cloud has a lifetime of approximately 
20-min during which these cloud particles will grow and perhaps evaporate again. 

Cloud physics research under zero-gravity (zero-g) or low-g conditions offers 
solutions to many of these problems. Under zero-g conditions, the experimenter can place 
a drop in a chamber and observe it through a microscope for long periods of time. The 
droplet can be frozen, thawed, and another drop propelled into it. The migration and 
collection ot particulates that may be near or around the drop can be observed. 
Numerous experiments (e.g., determining whether a freezing drop splinters and/or 
acquires a charge) can be done in this unique environment that cannot be done on earth. 
Other important experimental areas include the diffusional growth of drops and ice 
crystals and studies ot the effect of temperature on an ice crystal’s type and form. 

In the study of cloud physics one is primarily interested in and working with 
particles of some form and size. The physical parameters of the particles of interest arc 
many and varied. With regard to physical size, the particles of concern to the cloud 
physicist fall into basically two categories: (1) the condensation nuclei, which range from 
0.1 pm to about 1 pm, and (2) the giant nuclei, which include all particles greater than 
1 pm. The total particle size spectrum of interest is presented in Figure 1 . Figure 2 
presents some of the experimental chambers currently under consideration for the Zero-g 
Atmospheric Cloud Physics Laboratory proposed for the Shuttle/Spacelab Program. 

This zero-g-type atmospheric laboratory' has been conceived to provide the cloud 
physics scientific community with a general-purpose laboratory for furthering their basic 
understanding of cloud microphysical processes. During the course of the Zero-g 
Atmospheric Cloud Physics Laboratory feasibility study [1,2,3], it was recognized that 
techniques to photograph and observe the events occurring during the experiments must 
be identified and that studies should be performed to select the best technique or 
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Figure 2. Proposed atmospheric cloud physics chambers. 
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combination of techniques to accomplish the data management. Since holography has the 
unique capability to record large amounts of data that cannot be obtained by other 
techniques, a study was undertaken to determine which holographic system would be 
best for use with the Zero-g Atmospheric Cloud Physics Laboratory program. This report 
presents the results of the studies accomplished thus far. 


Specific Objectives of the Study 

The objective -of this study was to investigate the feasibility of applying 
holographic techniques to the real-time determination of particle parameters in a zero-g 
environment. Such feasibility would be judged, of course, in a 1-g laboratory 
environment, with the primary constraints of particle resolution and velocity being 
imposed on the proposed holographic technique. The particle parameters to be 
determined holographically were 

1. Particle Size 

2. Particle Size Distribution 

3. Particle Velocity of Motion 

4. Particle Velocity Distribution 

5. Particle Growth 

6. Particle Growth Rate 

7. Particle Trajectories 

8. Particle Deformation. 

If the environment or particle field of interest is static or dynamic and involves a 
reasonable volume, holography provides a tailor-made technique for the real-time 
detection, observation, and measurement of various physical parameters of the individual 
particles within the volume. Consequently, holography can provide much of the necessary 
information about the total content of the particle field. That holography is a tailor-made 
technique for such particle field investigations is evidenced by the fact that it allows the 
entire particle field volume of interest to be recorded on one exposure. Upon 
reconstruction of the hologram, one may inspect the recorded volume of interest for the 
investigation of the physical parameters of concern at any later, more convenient time. 
For example, consider that the astronaut experimenter has recorded a hologram during an 
experiment while in orbit. Upon reconstruction of this hologram in an earth laboratory, 
all observers would see exactly and precisely what the astronaut saw and recorded during 
the experiment. 
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Today there are numerous applications of holography to the problem of real-time 
particle detection. Some examples are fog and rain evaluation, wind-tunnel studies, 
rocket-exhaust analysis, pollution monitoring, and Skylab waste-tank analysis. Most of 
these previous applications employed the technique of in-line holography, which will be 
discussed in the following section. 


DESCRIPTION OF THE GENERAL HOLOGRAPHIC TECHNIQUE 


Introduction 

In general, if information on particles and particle parameters is desired from any 
given environment or particle field, one must resort to experimental techniques to 
determine particle parameter information. These environments or particle fields of interest 
may be anything from ordinary smokestacks, which produce common air pollutants, to 
the analysis of rocket exhaust, to wind-tunnel studies, to fog and rain evaluation, or to 
the present environment of particle fields within a controlled cloud physics chamber. 

As mentioned earlier, one of the first holographic techniques applied to such 
particle field investigations was the technique of in-line holography. With the aid of 
Figure 3, this technique may be described in the following way. A plane wave-front beam 
is incident on and passes through some target volume ot interest. The presence of 
individual particles distributed throughout this volume essentially causes the plane wave 
to become scattered and diffracted at many points so as to produce a spherical wave as 
well as a plane wave. The plane wave is, of course, that part of the incident wave which 
passes through the target field region undisturbed. This undisturbed plane wave then 
constitutes the holographic reference beam. The spherical wave produced by the 
scattering and diffraction around the individual particles constitutes the holographic 
object beam. These two waves interfere at the film plane, and the interference pattern so 
recorded at the film plane provides the means for reconstructing the real image of the 
target field region. An inspection of this real image provides quantitative information 
about the individual particles that make up the target field. The basic problem with the 
in-line system is that for high-density particle fields, the reference beam is excessively 
attenuated, and the observer or detector must look directly back into the laser source on 
reconstruction. 

The optical holographic system to be considered for application to the Zero-g 
Atmospheric Cloud Physics Laboratory program will be a modified acute sideband 
technique. Some basic descriptions of sideband holographic techniques, both for the 
stationary and moving object fields, will now be presented. These basic concepts are 
provided so that this study of an optical holographic technique can stand alone. 


5 




General Description of Sideband Holography for a Stationary Scene 

Let the plane of the hologram lie in the x-y plane of the coordinate system, as 
given in Figure 4. The amplitude in this plane can be described by the radiation from the 
scene, 


E 


s 


b (x,y) exp 


i[0(x,y) 



( 1 ) 


and by that of the reference beam, 


E 


r 


a(x,y) exp |i[a(x,y) 



( 2 ) 


The amplitude functions, b (x,y) and a (x,y), and the phase functions, a(x,y) and 
0(x,y), are considered to be real functions. This, then, is a scalar treatment that does not 
take account of polarization effects. However, it is always possible to split the scene wave 
into two components, one with the electric vector in the plane containing the electric 
vector and the direction of propagation of the reference ray and the other with the 
electric vector perpendicular to this plane. Only the parallel component contributes to 
the interference pattern on the hologram and, hence, to the reconstruction. Furthermore, 
a(x,y) is a slowly variable function (a constant for a plane reference wave but, for a 
spherical wave, it varies as the reciprocal of distance from the reference source to 
different points on the hologram). The reference wave may also be a distorted spherical 
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Figure 4. Coordinate system. 


wave, as formed by a lens with geometrical aberrations. The amplitude function of the 
scene (object) wave, b(x,y), is materially smaller than a(x,y) since the ratio of intensity 
of the reference beam to that of the scene beam is usually required to be from 3 to 1 to 
10 to 1. The exposure on the hologram plate is given by the product of the intensity at 
the plate times the time of exposure, or pulse length. 

In general, from the definition of Poynting’s vector in the rationalized 
meter-kilogram-second system 


7 


S = E x H 


one can see that its time average is 


<S> = (ExH) = 1/2 \fTfix EE* n 


( 3 ) 


where E* is the complex conjugate of E, and n is a unit vector, mutually perpendicular 
to i and H . Since the intensity (i.e., irradiance) of electromagnetic radiation is the 
energy crossing normal to a unit area per unit time, 


I = <S> = 1/2 y/e/fi EE* = 1/2 ZEE* = m EE* 


(4) 


where Z = \Zm/7 is the impedance of the medium, and m is a constant for simplification. 
For a hologram, 


I = m EE* = m (E f + E g ) (E* r + E* s ) 


(5) 


and from equations (1) and (2), 


I i[j3 (x,y) - cat] i[a (x,y) - wt] | 

I = m \b (x,y) e + a (x,y) e I 


, . . -i[^ (x,y) - cot] -i[a (x,y) - wt]| 

b (x,y) e + a (x,y) e i 


( 6 ) 


or 


1 = m ( b 2 (x,y) + a 2 (x,y) + a(x,y)b(x 


>y) [e* 


[(J (x,y) - « (x,y)] 


+ e 


-i[)3 (x 


,y) - a (x,y)]^ 


(7) 


where equation (7), for the intensity recorded at the plate, is generally called the 
recording equation. 
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Now the substitution of equation (7) into the equation for the exposure produces 


T T 

£ = f I dt = m / b 2 (x,y) + a 2 (x,y) + 2a (x,y) b (x,y) cos 03 - a) dt 
0 0 


or 


It = niT[b 2 (x,y) + a 2 (x,y) + 2a(x,y) b (x,y) cos (/3 - a)] (8) 

but r is just some constant exposure time, r 0 , for a given hologram; therefore, the 
exposure may be simply written as 

It = m t 0 [a 2 + b 2 + 2ab cos (j3 - a)] 


where a, b, a and j3 are coordinate dependent, or 


It = m t 0 


+ b 2 + ab 




i((3 - a) 


+ e 


-i (0 - «) 


(9) 


It is well known in photography that the Heurter-Driffxeld characteristic curve as 
plotted on the semilog plate is linear for most of its length and can be represented in this 
linear region by 


It 

D - y log — 
g 


where D is the density of the image on the film, y is the contrast, and g is the inertia 
of the emulsion. Since 


D = - log T , 

T = g 7 (Ir)“ 7 (10) 
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where T is the transmission coefficient. The amplitude transmission coefficient, T a , is 
the square root of the transmission; consequently, 


r yl 2 n - y I 2 _ • lc r yl2 

r a = g (Ir 0 ) - K1 


( 11 ) 


7/2 

where k = (g/r 0 ) is a constant and 


T a = k 


| f iff? - a) -i(f 3 - a)l 1 1,2 

I a 2 + b 2 + ab |^e + e J I 


( 12 ) 


Factoring, 


-7/2 ab i ((3 - 

T a = (k) (a 2 + b 2 ) 1 + 2 — e 

a a* + b 2 


a) — i (/3 - a)"| 

+ e j 


-7/2 


. ( 13 ) 


Binomial expansion may be used to obtain 


-7/2 I 7 ab i(0 - a) -i(0 - 

T, = « a = + b») ,-I_ e 


i ((3 - a) -i(/J - a) 
e + e 


+ higher-order terms 


( 14 ) 


However, by careful development of the exposed plate, one obtains a positive 
transparency with an overall gamma of (7 = -2), which yields 


ab 


Ta • W + »’> 1 ^ 77*. 


i(0 - «) -i(/J - 

e + e 


oj 


05 ) 


and higher-order terms are omitted because they are negligibly small compared with one. 
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Consider the usual case where the reference beam during reconstruction is 
essentially identical with that during recording. Let this reference wave be incident on the 
developed hologram; just beyond the hologram the wave amplitude is 


E t - T a E r 


(16) 


Then, from equations (2) and (15), 

i(a - to t) ab I i(0 - co t) -i(0 - 2a + co t) 


i t = Jka (a 2 + b 2 )J je 


a 2 + b 2 


i( 

e 


+ e 


■II- 

(17) 


The first term, 


ka (a 2 + b 2 ) e 


i(a - cot) 


(18) 


is then simply the reference beam attenuated by a term proportional to (a 2 + b 2 ); 
therefore, the first term is 


k(a 2 + b 2 ) E r 


(19a) 


The second term, 


[ ka(a!+b i[t?TF-) 


becomes 


i(/3 - co t) 

ka 2 (be = k a 2 E<. 


(19b) 


11 



which corresponds to the original wave, equation (1) multiplied by the same shading 
factor as for the reference beam. One may infer from the form of the second term that 
upon viewing through the hologram, illuminated by the reference beam, one should see 
this virtual image of the object in the precise position previously occupied by the object 
with reference to the original position of the hologram. If the hologram is in its same 
position after reconstruction, the image of the object will be precisely superimposed on 
the actual object. However, the brightness of the image will be slightly reduced. 

The third term, 


k b e ' iW ' 20 + “ l) 


i2 (a - (3) + i(0 - gj t) 

k a- b e 


becomes 


k a 2 


i2(a - (?) c 
e c s 


(19c) 


which corresponds to the real image of the object and is located on the opposite side of 
the hologram plate from the virtual image. This real image is not aberration free and is, 
therefore, termed pseud oscopic. 


Theoretical Description of Time-Dependent Holography 

The previous section was useful to show, mathematically, how a hologram is able 
to reconstruct a three-dimensional image of a real object by the reconstruction of the 
various wave fronts. We now intend to show that if the real object moves during the 
exposure of the hologram, the wave fronts become time dependent. For example, in the 
present problem of the zero-g cloud physics chambers, the particle field of the 
environment will be in constant motion. The effect of this object motion on the 
reconstructed wave front, i.e., the holographic image, must be determined. This 
information will allow determination of whether it is feasible to expect to record 
holographically such moving particle fields. 

From the previous section, the basic exposure equation may be written as: 


£(P) = m 



E 2 


+ 2 


E r E s cos k 


[s(P) - r(P)l > dt 


( 20 ) 
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It must now be determined whether the holographic system is to become time dependent 
or time independent. A time-dependent system will be considered first. This time 
dependence is permitted by having the object (e.g., a particle) moving during the 
exposure. Then during the hologram exposure, s(P), which is the object beam path length 
in equation (20), becomes 


s(P) = s p (t) 


( 21 ) 


and equation (20) for the exposure becomes 


&(P) = m t < 


F 2 + F 2 
c r c s 


+ 


2E,E 0 
r s 


r/2 

/ cos k [s (t) 
-r/2 



( 22 ) 


Therefore, the time dependence arrives as a result of the time-dependent path 
length change of the object beam caused by the motion of the object or particle. To 
proceed further, one must be able to completely specify the path length change s(t). For 
this discussion, consider the object or particle to have linear motion over some distance s 
during the time of exposure. Then, in general, 


ds 

dt 


and integrating, 


(23) 


s(t) = So + vt 


(24) 


The diagram for this very special case of the time dependence of s(t), along the line of 
the original s , is given in the Appendix, where the more rigoroui vector description of 
this time dependency of motion is presented. The purpose here is only to demonstrate 
the dependence of motion holography on the system geometry. Substituting equation 
(24) into equation (22) produces 


£(P) 


m t 


K C + 


2E r E* 


r/2 

j cos (kvt + 0) 
- r/2 


(25) 
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where = E f 2 + E s 2 and <i> = k (s 0 - r) =£ f(t). Performing the indicated integration of 
equation (25), 


£(P) = mr|K^ + 2E r E g sine kv r/2 cos <pj 


(26) 


where sine kv r/2 = [sin kv(r/2)] /[ kv(r/2)]. 

Equation (26) becomes the exposure equation for the time-dependent case of 
linear motion of a particle. Recall that equation (20) is the exposure equation for the 
time-independent case of the stationary particle. For this time-independent case, simply 
let 


[s(P) - r(P)] = 4> 


and integrate over dt to obtain 

&(P) = mr^K c + 2E r E s cos (27) 

as the final exposure relation for the time-independent case. 

If the particle velocity were zero, then the exposure expression of equation (26) 
should simply reduce to that of equation (27) for the stationary case. Setting v = 0 in 
equation (26), one need only treat 


sine kv r/2 = lim sine kv r/2 = 1 
v->0 


(28) 


and, indeed, equation (26) does reduce to equation (27) as expected. 

Therefore, the time dependence of holography is determined by the path length 
change of the object beam during the exposure. This change of path length is then 
reflected as an amplitude modulation of the phase term, cos <j> , in the exposure relation. 

It must be pointed out here that this change in path length of the object beam is 
a function of the direction of the motion vector for the object or particle with respect to 
the vector for the path length of the object beam. Because the path length change of the 
object beam results in a relative phase change between object and reference beam, the 
following section is presented. 
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Dependence of Relative Phase Shift on the Geometry of the 
Holographic Arrangement 

It has been shown that the time dependence of motion holography results from 
the path length change which occurs in the object beam because of the total motion of 
the object or particle during the exposure. Because holography is an interference 
phenomenon, the total change in path length, Ad , of the object beam during the 
exposure must be less than A/2 (i.e., kX/2 = jr). Because of this severe restriction, the 
geometry of the optical arrangement used becomes of utmost importance to the 
successful recording of a moving object or particle. Attention, therefore, is directed to 
two limiting cases of sideband geometry. 

Consider the holographic arrangement of Figure 5. This constitutes a very 
desirable arrangement for a stationary object because of the high energy_ return from the 
object to the film plate. For this case, the radiation propagation vector k is both parallel 
and antiparallel to the direction of the intended object motion, Ax. Consider that during 
the exposure the object moves from position x 0 to a new position x 1 ; we call this object 
translation Ax. Then the total optical path length change is 


Ad = 2 Ax 


(29) 


yet, using the limiting requirement for the allowed change in optical path length of a 
holographic arrangement, 


Ad < X/2 


(30) 


or 


Ax < A/4 


(31) 


and the object cannot travel a total distance greater than A/4 if a hologram is to be 
successfully recorded. This was an obvious result when the object motion, Ax, was 
allowed to be along the direction of the radiation propagation vector, k . This geometry 
allows the minimum amount of object motion during the exposure. 

A contrasting case is given in Figure 6. The basic difference in geometry between 
this and the previous arrangement is the rotation of one mirror. Everything is as before 
except that now the direction of motion, Ax , of the object is perpendicular to the 
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Figure 5. Maximum sensitivity sideband holographic arrangement. 


propagation vector, k . The result is, of course, that now the magnitude of Ax may be as 
large as desired without any change occurring in the object beam path length. For this 
geometrical arrangement, 


Ad s 0 


(32) 


Of course, if one expects to resolve the object in the hologram with such an arrangement, 
the allowed Ax is not totally unlimited. As a rule of thumb, it may not travel further 
than one-tenth of its length during the exposure because of resolution requirements. But 
this is an unusually large distance compared with X/2 . The price paid for this large 
motion is that no resolution of front-surface detail has been recorded. The image simply 
appears as a black silhouette against a lighted background. 

A unique geometrical arrangement has been conceived which provides the 
advantages of both of the previous limiting cases and apparently does so without the 
concomitant disadvantages of these systems [4—10]. This new system allows the 
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resolution of front-surface detail to be recorded from a moving object, gives maximum 
utilization of the available coherence length, and allows sufficiently large object motion, 
Ax , during the exposure to capture moving objects with velocities of extreme magnitude. 

Consider the ellipse of Figure 7. The three paths shown are all of constant length 
and equal to twice the value of the semimajor axis. This fact of constancy proves to be 
quite beneficial for the geometry. 

Suppose one could cause the object or particle of interest to move along the 
surface of such an ellipse. The laser source could then be positioned at one focus and the 
film recorder at the other. In this way, regardless of the total motion during the 
hologram exposure, the object beam path length would always be constant and equal to 
2a. Further, one could always front illuminate the object or particle. From a practical 
point of view, such an elliptical motion trajectory is not very useful. The situation 
previously described can be approximated, however, by allowing the' motion vector for 
the object to be tangent to the surface of the ellipse. If the motion is considered to be 
linear, then the linear motion vector must be parallel to the semimajor axis of the ellipse 
and perpendicular to the semiminor axis. In this way a straight line is used to 
approximate the surface of the ellipse. It turns out that the approximation is quite good 
over a small section of the ellipse surface [3] . This statement is verified by the ability to 
reconstruct a hologram of an object in constant motion during exposure [4] . This system 
has been used to successfully record a hologram of an object having total motion, Ax , 
equal to 1 0 000 wavelengths of the source radiation while allowing less than one-half 
wavelength changes in the object beam path length [9] . 


Some Technical Aspects of Image-Recorder Resolution and 
Effects of Recorder Size 

The point of view taken for the discussions here will follow that of Collier, et al. 
[11]. Consider the constant subject-reference phase difference of holography caused by 
the finite extension of the subject. The bright fringes of the interference pattern occur 
when 


AC = nXi 


(33) 


where n is an integer and AC is the constant path difference of point references 
between subject and reference beams. 

From the basic consideration of the concentric in-line system of Reference 1 1 , 
equation 3.5, 
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AC 


( 



where the numerical subscripts 1 and 2 refer to the subject plane and hologram plane, 
respectively, and the subscript r refers to some arbitrary reference plane in the reference 
beam, a distance Z r from the hologram plane. 

Substitution of equation (33) into equation (34) yields the general expression for 
a circle whose center has the coordinates 


Z 1 x r - Z r x 1 
Zj - z r 


and 


y 2 


z i y r - z r yi 
Z I - z,. 


(35) 


and whose radius P is given by 


p2 = 


Z 1 x r" z r X 1 
Z, - Z_ 


z i y r - z rvi 


Z 1 " Z r 



Z,, Zl Z| z r a 

Zj - z r 


(36) 


An acute sideband holographic system will be discussed later. Provided here is a 
discussion of just what the effect of such an acute angle means with regard to the 
interference pattern recorded. 

Consider Figure 8, which shows an off-axis intensity pattern formed by the 
interference of an axial plane reference wave (x r =y r =0, Z r = °°), with a spherical 

subject (i.e., object) wave diverging from an off-axis point (x, y = 0, Z). The center of the 
set of circular fringes, whose radii correspond to integral values of n in equation (36), is 
given by equation (35) as x 2 ' = x l an d y 2 ' = 0- The result is a zone plate pattern 

centered at the foot of the perpendicular dropped from P to the plane of the hologram. 
From Figure 8, if the photographic plate is centered at zero so that it can record only 
off-center portions of the interference pattern, a sideband hologram is obtained. If 0 j of 

this Figure is zero, the in-line situation results. For 6 j equal to a small acute angle, the 

acute sideband holographic system is formed. The fringe frequency £' in 
the x 2 ' direction can be found by differentiation of AE/X of equations (33) and (34) 

combined, under the conditions x r = y r = y j = 0 and Z r = <*> . The result is 


Z, X 


1 A 1 


(xi - x 2 ') 


(37) 
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Figure 8. Holographic zone plate produced by off-axis object P. 


For the in-line arrangement, the fringe frequency from Reference 1 1 is given as 


«i' 


~ x 2~ 

z i M 


( 38 ) 


As before, assuming the hologram plate to be centered at zero, a comparison of 
the fringe frequencies for the two arrangements may be accomplished. At the hologram 
center (X 2 ' = 0), the in-line fringe frequency is zero, while the off-axis fringe frequency is 
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xj/ZjXj . As one proceeds outward in the negative X 2 ' direction (Fig. 8), the frequency 
of each fringe system increases linearly with X 2 * and this frequency difference is 

maintained. At the edge of the hologram are generated the highest frequencies to be 
recorded. If the photosensitive medium on the hologram plate is to record the fringes in 
the off-axis case, it must have a resolution capability Xj/ZjAj in addition to that 

required for the in-line case. An examination of equation (34) reveals that if an otf-axis 
plane wave were considered (x r # 0, x r '/Z r = tan 0 r « 0 r ), the frequency difference 
would be 

1 (Hi . — 

A 1 \ Z 1 Z r 

where 0j is the mean angle to the axis made by the subject wave, i.e., the angle to the 

z-axis made by a ray passing from P to the center of the hologram at 0. Thus, the mean 
angle between subject and reference beam provides the difference in the maximum fringe 
frequency generated in an off-axis hologram compared with the in-line hologram. 

In the practical case of an extended subject, either the width of the subject or 
that of the recording plate can cause the fringe frequency i to exceed the plate 
resolution R . When the plate is small compared with the subject, then the 
term Xj/Zj Xj or Q ] -0 r /\ j is dominant. Point sources at the extreme dimension of the 

subject produce the maximum fringe frequency at the plate. If g'>R tor the extreme 
portions of the subject, then such portions are not recorded. On the other hand, when 
the plate is much larger than the subject, the dominant term in equation (37) might be 
the first. Beyond some value of x 2 ' all subject points produce zone plate fringes whose 

frequency %' > R . That value of x 2 ' defines the practical extent of the hologram record. 

Limited recorder resolution will now be discussed because ot its importance to 
the problem at hand. Consider the x-z plane of Figure 9, where a subject, a reference 
point source R(x r , 0, -Z r ), and a recording plate are indicated. The plane normal to the 

axis and containing the reference point source is separated from the recording plate by 
the axial distance Z r . As indicated in Figure 9, a ray from the reference source to an 

arbitrary point Q(x 2 , 0, 0) on the plate makes an angle & r with the normal to the plate 
(the z-axis) and is, therefore, characterized by the spatial frequency £ r Suppose that a 

small bundle of reference rays interferes in the vicinity of Q with a corresponding bundle 
of rays from the subject making an angle 0 q with the z-axis and characterized by a 

spatial frequency £q. The interference pattern intensity which must, be recorded is 


\ _ e _l “ e r 
/ ~ 
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Figure 9. Recorder resolution geometry. 


I - £exp ( 2 n i £q x 2 ) + exp (2w i | r x 2 )J 

• £ exp (-2ir i | Q x 2 ) + exp (~2n i £ r x 2 )J 


= 2 + 


2 




(40) 


where unit amplitude waves have been assumed and the fringe frequency is (|q - £ r ). For 
small angles, 


sin e q - sin 6 r 

^0 " ^r “ s 



( 41 ) 
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Suppose now that the recording medium is perfectly capable of recording all 
fringe frequencies below a cutoff frequency, £ c , but is totally incapable of recording 

fringe frequencies exceeding | c . For a fixed value of | r and a subject which scatters 
light over a broad band of spatial frequencies £q , there will be limits to the subject 

spatial frequencies that can be stored in the hologram; that is, it rays from the subject 
are to be stored in the hologram, their angles with the z-axis must lie within certain 
limits. In Figure 9, one such limiting ray passing from the subject to an arbitrary point Q 
on the hologram plane intersects the plane of the reference source at the point x max . 

We call this ray a marginal ray, and its spatial frequency £q max satisfies 


£o,max ^r ^c 


(42) 


If the angle the marginal ray makes with the z-axis is 0 Omax > then from equation (41), 


^0,max ^r 
X 



(43) 


or 


0 O,max " 6 r 58 ^c 


(44) 


From Figure 9, the point x max can be obtained and is given by 


x max ~ x 2 

zT“ 


tan (e r - e c ) d T + d c 



£ c X 


or 


x max 


x r + Z r X £ c 


(45) 
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Likewise, for 


x min 


Y Y 

min 


z r A ?c 


( 46 ) 


Note that the marginal rays intersect the plane of the reference source at a coordinate 
independent of the point Q in the recording plane. Note also that if the recording 
material will not record spatial frequencies greater than j| c )o (the absolute value of the 

cutoff frequency), the resolution in the image is limited accordingly. This will be of 
utmost importance in choosing the best film recorder. 


ANALYTICAL DESCRIPTION OF THE ELLIPTICAL HOLOGRAPHIC 

ARRANGEMENT 


The system to be presently discussed provides the basic system which, when 
modified, will furnish the technique for use as a real-time measurement system for the 
Zero-g Atmospheric Cloud Physics Laboratory investigations of moving particle fields. 


Basic System Description 

The elliptical technique is so called because of the relative location of. the various 
optical components with respect to each other. Depending on the magnitude of the 
expected object motion, a very specific set of locations for the optical components is 
chosen (i.e., a very specific elliptical arrangement). Further, depending on the direction of 
the object motion, a very precise orientation of this elliptical arrangement is chosen with 
respect to this motion vector. 

The specific orientation of the holographic system just referred to is based on the 
use of the specific elliptical holographic arrangement oriented with its major axis parallel 
to the motion vector defined by the moving object. This motion vector must be made 
tangent to the mathematical surface of the specific ellipse at some point Q. A typical 
configuration of this holographic system positioned in the preferred orientation is shown 
in Figure 10. 

This unique system is defined by the following conditions: a thin film beam 
splitter (B/S) centered at the focus, fj, of the chosen ellipse; a film recorder centered at 

the other focus, f 2 ; and the major axis of this ellipse, defined by xx' in Figure 10, being 

parallel to the tangent line, PP' , which is identified as the motion vector of the 
high-speed particle or object. Laser radiation is incident on the thin film beam splitter, 
centered at focus fj . The transmitted beam from here is made incident on the object or 

particle which is translating linearly along the tangent line, PP' , and is momentarily at 
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the point, Q, on the perpendicular bisector of xx' at a distance, b, from the major axis, 
where b is the value of the semiminor axis for the chosen ellipse. This object beam is 
then made incident on the film positioned at the focus, f 2 . The beam reflected from the 

beam splitter at fi constitutes the reference beam for the system and, after a reflection 
from a mirror, is made incident on the film at focus f 2 . These two beams interfere at f 2 , 

and the hologram is recorded. Because of the judicious choice of paths for these two 
beams, their path lengths are matched, and maximum utilization of available coherence 
length is accomplished. 



Figure 10. Typical elliptical holographic configuration. 


Linear Motion of Object in Terms of Parameters of the Ellipse 

Consider Figure 1 la; the general equation of such an ellipse is given by 


b 2 x 2 + a 2 y 2 = a 2 b 2 


( 47 ) 


The line segment, PP' , is considered to be tangent to this ellipse at the point Q, which 
lies on the perpendicular bisector of xx' . This line, PP' , in Figure 11a is identical to the 
tangent line, PP' , of Figure 10. It is the so-called line of motion of the high-speed object, 
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Figure 1 1 . Family of ellipses. 

is parallel to the major axis, xx' , of the ellipse, and may be considered perfectly straight. 
The projectile travels along PP' of Figure 1 la and reaches point Q at some time tQ . The 

radiation incident at this point at tQ will be reflected to the film, which is positioned at 
focus f 2 . At this particular moment the mathematical ellipse passes through point Q, 

with a beam splitter at fj and a film at f 2 , and results in the situation depicted earlier 
in Figure 10. 
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As the object moves some incremental distance, Ax , along PP' , past point Q, it 
moves off this initial ellipse, but it can be considered to move immediately onto another 
ellipse just slightly larger than the initial one. If the elliptical constant of the initial 
ellipse was 2a, then the elliptical constant of this new ellipse will be 2(a + Aa). The 
radiation reflected from this moving object will then be incident on the film at f 2 and 

will interfere there with the reference beam, as long as 2 Aa is less than A/ 2 . 

In Figure 1 lb a family of such ellipses is constructed, each successive ellipse being 
intercepted by the line segment, PP' , as one moves from Q to the right along PP' 
parallel to the x-axis of the coordinate system. The separation of the foci must remain 
constant and equal to 2d for the entire family of curves. Figure 11c is a convenient 
enlargement of the first quadrant of Figure lib. The points of interception of PP' with 
each successive member of the family of ellipses are given by P j , P 2 , P 3 J etc., 

respectively. As PP' is traversed to the right, the original ellipse can be considered to 
grow successively to the next larger member of its family, while the foci separation 
distance 2 d remains constant. 

If one considers that the ellipse is to enlarge during the time t, then expanding 
to a larger ellipse, equation (47) becomes 


b 2 x Ax + b Abx 2 + a 2 y Ay + a A ay 2 = a 2 b Ab + b 2 aAa (48) 


Since the direction of object motion is parallel to the x-axis, Ay = 0 and 


b 2 x Ax = a A a (b 2 - y 2 ) + b Ab (a 2 - x 2 ) 


(49) 


However, from Figure 10a, 

a 2 - b 2 = d 2 , (50) 

where d is a constant; therefore, 
a Aa = b Ab , 
and equation (49) becomes 

b 2 x Ax = a Aa[a 2 + b 2 - (x 2 + y 2 )] . (51) 
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From the basic equation for the ellipse, equation (47), one may easily obtain 


x 2 + y 2 = b 2 + x 2 [ 1 - (b 2 /a 2 )] 


( 52 ) 


Substituting equation (52) into equation (51) gives 


x Ax = (Aa/ab 2 ) (a 4 - d 2 x 2 ) 


(53) 


For any ellipse, 2a = L, where L is a constant and at present is the optical path length of 
the front illumination signal arm. When the ellipse expands as a result of object travel 
along PP', then 2a= L changes by 2Aa = AL. Therefore, since Aa = AL/2, equation (52) 
becomes 


x Ax = (AL/2) (1/ab 2 ) (a 4 -d 2 x 2 ) 


(54) 


Now since point Q has been taken as the reference point for x (i.e., x is zero when the 
projectile is at point Q), as measurement of the projectile motion starts from 
point Q and traverses to some point P, x = Ax and equation (54) becomes 


(Ax) 2 = (A L/2) (1/ab 2 ) [a 4 - d 2 (Ax) 2 1 


(55) 


or 


(Ax ) 2 


AL a 3 /2b 2 
1 - (ALd 2 /2ab 2 ) 


(56) 


Making the assumption that the magnitude of d 2 is not drastically different 
from b 2 and AL/2 « a, 


AL d 2 /2ab 2 « 1 


and equation (56) becomes 


(Ax ) 2 = (AL/2) (a 3 /b 2 ) 


(57) 



where AL = 2Aa is the variation in the original elliptic constant 2a resulting from the 
travel of Ax of the object along PP'. At a given velocity V for a time t, 


Ax = Vt 


(58) 


where t is the pulse length and also the exposure time. 

Equation (57) is, therefore, an expression which relates the distance Ax traveled 
by the object along PP' to the total change in tile elliptic constant, AL = 2Aa; i.e., AL is 
the change in the optical path length of the front-surface illumination arm of the 
holographic system. Substituting equation (58) into equation (57) and solving for the 
object velocity, 

V = (AL/2)’^ a 3/ Vbt . (59) 


This relation may be used to determine the permissible object velocities allowed by the 
specific configuration having a set of elliptic parameters and a specified tolerance AL. For 
illustration, one may arbitrarily set AL equal to X/8[X = 694.3 nm (6943 A)] and let the 
distance of separation between the center of the first beam splitter, at f j , and the center 

of the photographic film, at f 2 > be a constant value 2d. By varying the semimajor axis a, 

which in turn varies the semiminor axis b, a set of permissible velocities can be obtained. 
This set of velocities is graphically shown in Figure 12, where the permissible velocity 
values have been used as the ordinate and the arbitrarily chosen values of the semimajor 
axis a as the abscissa. Each separate curve corresponds to a specific value of d, and the 
elliptic parameters are related by a 2 = d 2 + b 2 . 

It may be interesting to note the following three items. First, for this fixed value 
of A L and each assigned value of d , the curve approaches the vertical line, a = d, 
asymptotically. This seems to indicate that the projectile velocity can be any high value 
without limit if a = d. Obviously this is not practical since at a = d, b = 0, and the object 
would have to pass directly through the beam splitter and film. However, picking the 
smallest practical value of b allows the highest possible velocity for a given value of d. 
As the assigned value of d increases (bounded by some practical value of d), the curve 
rises and thereby raises the allowed value of velocity. However, because of the steepness 
of the curve, this region may be somewhat unstable with respect to changes in a or b . 

Second, as the assigned value of d decreases, the respective curve lowers. The 
lower bound for these curves occurs at d equal to zero, i.e., the mathematical ellipse 
becomes a circle. This is again impractical because the beam splitter would be located at 
the photographic plate. 
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Figure 12. Graph of permissible velocities vs parameters of ellipses. 


Third, the differentiation of equation (59) shows that each curve has a minimum 
at the value ol a = \/3cT. Substitution of this back into equation (59) produces 


V min = (3 3/ */2) (dAL/t 2 )^ 


(60) 


This V min is tlle minimum permissible value of the velocity for each specific value of d. 

Because each curve has a zero slope at this point, the region about this point is the most 
stable region with regard to possible changes in the values of the elliptic 
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parameters a or b . (Changes in a, and therefore b, will occur as the mathematical 
ellipse enlarges because of object travel Ax; further changes in b, and therefore a, might 
occur because of the projectile varying slightly off path as it travels along the line of 
motion.) Because of this stability, this region probably constitutes the best choice for an 
actual experiment. 

In view of the constant consideration given to the direction necessary for the 
velocity vector for this system and the stated fact that this velocity vector must be 
parallel to the major axis of the chosen ellipse, it remains now to discuss the total depth 
of volume allowed by this system. 

An experimental and theoretical discussion of the allowed volume is presented in 
References 5, 6, and 11, where it is shown that the allowed volume is a function of the 
total displacement of the object during the exposure, Ax, and the parameters a, b, 
and d for the chosen ellipse; i.e., 


Vol = f (a, b, d, Ax) 


It obviously will be a function of the total volume originally illuminated since the 
enlarging or spreading of the illumination beam allows the system to function over a 
family of chosen ellipses instead of just one. This is explained by the fact that each 
separate angle formed by an enlarged or spreadout beam selects a separate member of the 
family of chosen ellipses. 


Consideration of the Primary Parameter of Motion Holography — 

Total Motion of the Object Field During Exposure 

The parameter of importance to motion holography is neither the velocity of the 
object, V, nor the exposure time, t, by itself. Rather, it is the constant product of these 
two parameters, which is the total distance. Ax, traveled by the object during the 
exposure time; i.e., 


Ax = Vr 


(61) 


This obvious statement has been reiterated for clarity, and Figure 13 is presented as a 
more useful engineering diagram for employment in developing the acute, sideband 
holography technique described later. This figure is only a modified version of Figure 12 
and may be used to determine the feasibility of obtaining a hologram of a particle 
traveling at 6.096 x 10 s cm/sec (20 000 ft/sec). 


31 



(uirt) XV 



Figure 13. Engineering diagram for determination of permissible object motion. 


From the experimental evidence of References 5 and 9, the value of total particle 
travel, Ax, during the exposure is allowed to be 700pm for a successful hologram. Using 
this value for Ax, Figure 13 indicates that an ellipse with a value of 50 cm for d will 
suffice. Then for any ellipse. 


a 2 = d 2 + b 2 


(62) 


and, depending on the geometry, a value of b as small as is practical is selected. 
Equation (62) provides the magnitude for the semimajor axis a. A specific ellipse has now 
been chosen by first obtaining an expected Ax which allows a successful hologram- 
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One must now use Ax and the expected particle velocity to obtain an exposure 
time or pulse length and then decide if this exposure time is realistic or feasible; if so, 
then the chosen holographic system is feasible. Consider that 


Ax = 700 x 10' 4 cm 


and 


v = 20 000 ft/sec = 6.096 x 10 s cm/sec 


then 


r 



1 1 4 x 1 0" 9 sec 


and 


t = 100 nsec 


(63) 


Therefore, this chosen system would be feasible with a ruby or pulsed argon laser as a 
source. 


BASIC DESIGN OF THE OPTICAL HOLOGRAPHIC TECHNIQUE PROPOSED 
FOR THE REAL-TIME MEASUREMENT OF THE PARTICLE 
FIELD GROWTH RATES 


The basic design to be used for the laboratory investigation of the various cloud 
physics chambers will be a modification of the basic elliptical configuration discussed 
previously. This modified basic design will constitute a laboratory-type system 
constrained only by the primary parameter of motion holography — the total motion of 
the particle field during the exposure and the resolution limit of the system. 

Consider Figure 14, which shows the modified basic elliptical technique (acute 
sideband). Based on the information of Figure 13, the parameters for this ellipse are 
given by 
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d = 50. 0 cm 


Figure 14. Modified basic elliptical technique. 


semimajor axis, a = 51.9 cm, 
semiminor axis, b = 13.9 cm, 


and 


(64) 


foci-origin separation, d = 50 cm 


Now from a combination of equations (58) and (59), an expression for the allowable 
total travel, Ax, is 



Now consider the quantity 

A L = 2 Aa , (66) 
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which is the total change in the optical path length of the front-surface illumination arm 
(object beam) caused by the total motion of the object during the exposute. This 
quantity, AL, must be less than X/2 for a successful hologram. The value 


AL = X/4 


(67) 


is chosen because this should be a good reconstruction of the moving particle. 

Substituting the values from equations (64) and (67) into equation (65) allows 
the calculation of the total expected travel, Ax, permitted by this system. 

■ -(± »V ;2 < 51 ' 9 cm) 3/2 (68) 

\ 4 2/ 13.9 cm 


and X for an argon laser is 5145 x 1CT S cm; therefore, 


Ax - 


f \ 1/2 3/2 

5145 x 1(T 8 cm \ (51.9 cm) 


8 


13.9 cm 


(25.359 x 10~ 4 cm 1/2 ) (373.896 cm 3/2 ) 
13.9 cm 


(69) 


Ax = 682.132 x 10' 4 cm 


Consider now a shutter speed of 1/60 sec commensurate with the lOfps of the Hulcher 
70 film transport. Then the useable pulse length is 


t = 17 msec 


r = 17 x 1 O' 3 sec 


(70) 
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Using 


£ - V 


(71) 


v = (682 x 10‘ 4 cm)/(17 x 10‘ 3 sec) 


v = 40 x 10' 1 cm/sec 


(72) 


v = 4 cm/sec 


The exposure time of 1/60 sec represents the longest exposure time even anticipated. 
Consider now the shortest exposure time which is readily available, 


t = 25 x 10" 9 sec. 


With this exposure time the allowable velocity becomes 


v > 31.6 96 x 10 9 cm/sec , 


v « 2 1 000 mph 


(73) 


Primary Constraints on Prototype System 

As mentioned, the primary constraints on the prototype system will be total 
motion during exposure and resolution limit. From equation (69) the total motion of 
682 pm is obviously more than adequate when recorded with the proper exposure time. 

The second constraint of resolution limit appears to be a more limiting one in 
view of the requirement of the smallest particle sizes to be observed in the zero-g 
atmospheric cloud physics chambers. The desired range on the particle diameters, Pq, of 
these chambers is 


0.01 pm < Pp) < 10 000 pm 


(74) 
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Since this range covers six orders of magnitude, no single system could be expected to 
monitor the entire range. The proposed technique, however, appears to be able to cover 
four orders of magnitude, 


1 Mm < Pj) < 10 000 pm 


(75) 


and still be able to reconstruct an image. Using basically the same system but reading the 
data in a different fashion, it is conceivable that the other two orders of magnitude could 
be obtained; i.e., 


0.01 pm < Pp < 1 pm 


(76) 


The details of tills will be covered in the next section. 


Analytical and Experimental Testing of Primary Constraints and 
Indication of System Response 

With regard to the primary constraint of total motion during the exposure of the 
hologram, there appears to be no problem as long as this motion is approximately 500 to 
700 pm. This total motion can always be controlled in the present situation of the cloud 
physics chambers by proper control of the pulse length or exposure time of the laser 
source. 


Several tests have already been performed in an effort to ascertain the total 
allowed motion. Figure 15 displays the experimental arrangement in which the free-fall 
motion of calibrated glass bead particles was used. Using the total allowed motion of 
700 pm and g = 980 cm/sec 2 , the exposure time was properly adjusted and the spherical, 
calibrated particles were successfully recorded. 

Using the same system as shown in Figure 15, the hopper and calibrated particles 
were replaced by a water chamber (Fig. 16). Air was constantly fed into this water tank, 
and the rising bubbles of varying size were recorded. Figure 17 displays the results from 
this experiment. The chamber and bubbles were recorded with a lens interposed between 
the chamber and the film for prerecorded magnification (Fig. 17a) and without lens 
magnification (Fig. 17b). 

With this system and premagnification, particles in the range of from 1 to 3 pm 
may be recorded and imaged. The relation between recording distance and particle size 
for this system was established by using various arrays of particles and a high-resolution 
test target. 
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Figure 1 5. Experimental arrangement for holography of free-fall objects. 



Figure 16. Experimental arrangement using water chamber. 


Figure 18 presents this graphical relation. Figure 19 displays a photograph of the 
virtual image of the high-resolution test target using premagnification. The resolved line in 
the lower right-hand corner is 251.2 lines/mm, which is roughly equivalent to the 
resolution of a 3-pm particle. 

These rough tests have been only for the basic response of the system. The 
intentions are to properly refine these techniques, apply both pre- and post-magnification, 
and automate the reading and detection techniques. Further intentions are to use the 
various cloud physics chambers for real-time detection of the existing particles. 
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a. With lens magnification. 



b. Without lens magnification. 

Figure 17. Experimental results using water chamber. 
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RECORDING DISTANCE (m) 



PARTICLE SIZE km) 


Figure 1 8. Engineering diagram for particle resolution vs recording distance. 



Figure 19. Photograph of high-resolution test chart recorded holographically. 
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This same system can be employed to detect the lower range of particles, i.e., 


0.01 jum < Pd < ljum 


if the data are handled in a different fashion. 

The image of such a particle is not reconstructed; rather, the hologram is treated 
as an interferogram. The use of this technique is called submicron holography and may 
be described as follows. The basic technique is to use the hologram as an interferogram. 
Then, from the intensity distribution of the diffraction pattern of an individual particle, 
one determines the radius, r, of the airy disc, This is related to the particle size by 


2a = 


1.22Z X 


where 2a is the diameter of the particle, and Z is the distance of separation of the 
particle from the film plane, i.e., object distance. The above equation is, of course, for a 
particle with a circular cross section. However, for a particle with an arbitrary cross 
section, one simply considers 2a to be the largest dimension of the particle. 

The limitation here is, of course, that there is no image with which to work. 
Obviously one can get size and size distribution trom this information, but any other 
information would be difficult because one could only monitor the motion of an array of 
airy discs caused by the presence of the particles. This technique must be investigated 
further. 


CONCLUSIONS AND RECOMMENDATIONS 


Four important processes in clouds must be better understood before deliberate 
weather modifications can occur. These are nucleation, growth, scavenging, and electrical 
charge separation. This necessary understanding can be greatly facilitated through the 
Zero-g Atmospheric Cloud Physics program. In the entire study of cloud physics, primary 
interest is centered on observing various parameters of particles having a wide range of 
forms and sizes. More importantly, the forms and sizes occur over and within some 
volume, i e., the various cloud chambers. Obviously, the simultaneous recording of a 
depth of particles is to be highly preferred over the recording of a limited planar cross 
section through this volume. This is the case of holography versus photography (Table 1). 
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TABLE 1. COMPARISON OF HOLOGRAPHIC AND PHOTOGRAPHIC RECORDINGS 


Characteristic 

Photographic 

Holographic 

Cost of Media and 
Replication 

Expensive silver halide 
media and wet processing 
required for archival 
master and hard copies. 

Ordinary vinyl tape suitable 
for a copy medium using 
hot embossing technique. 

Reduction Ratio 

Limited by scratches, 
dust, lens depth of 
focus (which is inversely 
proportional to the lens 
f/no.), and registration 
requirements of the 
system. 

Relatively immune from 
scratches; a large effective 
depth of focus, and 
relatively tolerant of 
registration errors. 

Brightness of Image for 
a Projection System 

Dependent on the 
luminance which the 
film can produce 
before distorting due to 
heat and proportional 
to the lens f/no. 

Dependent on source 
brightness and efficiency 
of hologram. (A phase 
medium does not absorb 
light, i.c., heat.) 

Optical Limitations for 
High Reduction and 
Display Resolution 

1. A large f/no. is 
required (to satisfy the 
diffraction limit 
requirements on the 
lens) which, for 
moderately priced lenses, 
is specified only for a 
limited field of view 
and a shallow depth of 
focus. 

2. High-cost lenses 
required to record and 
reproduce. 

1. Large depth of focus 
inherent in hologram. 

2. Requires lower-quality 
lens, especially if matched 
lenses are used to record 
and reproduce. 

3. Field of view limited by 
resolution of recording 
medium. 

Index Information 
Registration 

Misregistration causes 
loss of information. 

Misregistration causes loss 
of brightness. 
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The particles of concern to the cloud physicist fall into two basic categories, the 
condensation nuclei, which range from 0.01 pm to about 1 gm, and the giant nuclei, 
which include all particles greater than 1 jam. It has been shown that the entire particle 
size spectrum is from 0.01 uni to 10 000 jam or 1 cm. 

It has also been demonstrated that holography provides an image of the total 
volume of particles with one exposure. Therefore, holography captures the entire field 
distribution of particles within the total chamber volume instantaneously with one 
exposure. However, its resolution is limited to the top two-thirds of the particle size 
spectrum. Photography, on the other hand, does not have this resolution limitation, but 
it is severely limited by its inability to record the total volume. It must record only a 
planar cross section through the chamber. 

To properly utilize the advantages of both of these recording techniques, 
holography should be applied to all cloud physics chambers in the system described 
earlier. This will provide volumetric data throughout each of the individual chambers with 
regard to particle growth, trajectory, growth rate, and the other pertinent parameters. 
Additionally, for information on the particles falling in the lower one-third of the size 
spectrum, photography should be utilized. The two recording techniques can easily be 
employed in a hybrid arrangement. This hybrid system could be fashioned similar to the 
following. Use the modified basic ellipse system discussed earlier for the holographic 
system. Record the hologram on holographic film using either a 70- or 35-mm film 
transport. To employ the photography, rotate a turret to properly place the lens system 
for the camera which will utilize the same film and film transport. The illumination for 
the photographic camera system would be due to the object beam of the holographic 
system alone, i.e., use only one of the two beams in the holographic system. To use the 
holographic camera again, simply rotate the turret again to remove the lens from the 
object beam, and the holographic system would still record on the same film transport. 

Another possibility would be to use a beam splitter to allow both holography and 
photography to record events in a given chamber simultaneously instead of sequentially, 
as indicated above. In any case, a hybrid technique such as this would utilize the 
advantages of both recording techniques without any of the concomitant disadvantages of 
either. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Ala. 35812, Oct. 30, 1973 
645-10-01-000 
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APPENDIX 


THEORETICAL DEVELOPMENT OF HOLOGRAPHY 


Effect of Linear Scene Motion During Hologram Exposure 

The effect of scene motion during the exposure of a hologram is a spatial 
modulation of the recorded fringe contrast. This, in turn, causes a spatial amplitude 
modulation of the reconstructed wave front that blurs the reconstructed image in the 
same way that it would blur a conventional photograph having the same exposure time. 
In addition, this modulation reduces the brightness of even that part of the image that is 
reconstructed. The amount of object motion that can be successfully tolerated during the 
hologram exposure certainly depends upon the geometry of the holographic configuration 
with respect to the velocity vector of the motion. 

Two different theoretical approaches to the time-dependent theory of motion 
holography will be presented which follow closely the methods presented by Neumann 
[12] . Both approaches serve a purpose. The first uses a coordinate system with its origin 
at the scene point of interest, which is a vector approach that allows a more intuitive 
feeling for the physics of the problem. It is primarily valid only for those cases where the 
allowed motion is small compared with the distance from the scene to the hologram. The 
second approach uses a coordinate system with the origin at the hologram. It is the most 
general derivation, and if a sufficient number of terms are taken in the infinite series that 
occurs, it will cover all geometries. It suffers, however, from the mathematical complexity 
that results from these expansions. 

Scene-Oriented Coordinates. Following Neumann, a basic holographic 
configuration is assumed (Fig. A-l). The radiation from the laser is split into two beams. 
As usual, these are considered to be the reference beam and the signal beam. Let the 
point R of Figure A-l serve as the reference position for the system and S be some point 
of the scene. Light scattered from point S is made incident at the point P of the 
hologram, where it now interferes with light from the reference point. 

Let the field at point P be given by 


or 


E(P) = E r (P) + E S (P) 


E S (P) 


i[wt - ks(P)) 


E(P) = E r (P)e 


i[wt - kr(P)] 


+ 


(A-l) 




where E r (P) is the contribution from the reference point and is considered to be plane 

polarized, E g (P) is the contribution at P because of the scene point, S, and is the 

component of the scene-field parallel to the reference polarization, and r and s are the 
corresponding path lengths from the laser to the point P. 

Now, the intensity is given by 
I = m |E| 2 


Therefore, 


I = m lE r 2 + E s 2 + 2E r E s cos k(r - s)j 


(A-2) 


If the point, S, is allowed to move with some velocity, then s becomes a function of 
time. It is assumed, however, that the field amplitude does not change appreciably over 
the time, r, of the exposure. Thus, the interference pattern remains approximately 
constant except for the motion of the fringes. 
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Then, as before, the exposure can be written as 


r/2 f 

£ (P) = m f I E f 2 + E s 2 + 2E r E s cos k[s(t) - r] 
-r/z i 


dt 


and 


£(P) = m r 



F 2 


2E r E s 


r/2 I 

y cos k[$(t) - r] > dt 
-r/2 


Let 


K C = E r 2 + E s 2 


and 


&(P) = mr 



2E r E s 


r/2 

/ 

-r/2 


cos k[s(t) - r] dt > 


(A-3) 


To proceed, the time dependence of s must be given. In general, 
s = Sq + f(t) where f(t) may be any function of time where the midpoint of the exposure 

is at t - 0 . 

In particular, if it is assumed that s varies linearly with time, the result is 


s = Sq + vt 


then equation (A-3) becomes 


£ (P) = in r 




2 


r/2 

J cos [kvt + k(SQ - r)] dt 
-r/2 


(A-4) 
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Let 4> = k(sg - r) be the time-invariant portion; then 


& (P) = m r 


K c + 


E r E s 

T 

2 


t/2 

J cos (kvt + 4>) dt 
-t/2 


(A-5) 


after integration: 


£ (P) = in r 



E r E s 


sin kv - 

2 



cos $ 


(A-6) 


Then, since sinx/x = sine x, one obtains for the exposure 


£(P) = mr[K(- + E r E § sine (kv r/2) cos $ ] 


(A-7) 


If V = 0, then no motion exists (i.e., the stationary case), and the sine function 
becomes unity. This is the most desired case for most holographic experiments, i.e., that 
the argument of the sine function be zero. 

It is observed that if V is constant over the hologram, the resulting fringe 
contrast of the exposure will be poor for large values of the sine argument, and the 
resulting fringe contrast will be zero whenever the argument is mr(n = ±l, 2 ...). Consider, 
for example, n = 1, where the first zero occurs; then, 


kvr 



(A-8) 


or 


VT = X 


(A-9) 
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Now, since vr is the change in the distance, s, during the exposure, the 
possibility exists that motions of the scene as small as a wavelength may totally destroy 
the fringes or, at best, result in a poor hologram. Fortunately, V is not the same for all 
points of the hologram, and the fringe contrast is spatially modulated rather than 
destroyed. 

This special case of the linear time variation of s has demonstrated the sine 
function modulation of the recorded fringes. The problem may now be discussed with 
respect to the velocity vector of the motion. Consider the diagram in Figure A-2, where 
R is the reference point and, as before, S is the scene point. Let L be a point on some 

reference phase front from the laser, which will be common to both beams. Let a^ be 
the unit vector along the direction of propagation of the incident laser radiation at the 
scene point, S; £L be a unit vector along the line from scene point, S, to hologram point 

r 

P; and be a unit vector indicating the velocity of the scene point, S. 


LASER 



The reference beam is constant in length and composed of two parts, LR + RP. 
The front illumination beam or scene beam is composed of the two parts, LS + SP. The 
rate of change of the path length, LS, from the laser to the scene point, S, is given by 


V 



V cos 7 


(A- 10) 
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The rate at which the path length, SP, from the scene point to the hologram point, P 
changes is given by 


V • * 


a„ = v cos 6 


(A-l 1 ) 


Therefore, the path length from the laser to the hologram point, P, through scene 
point, S, is increasing at the rate 


ds _ a ^ A 

— - (V • a k - V • a p ) 


(A- 12) 


/ ds = 
s 0 


r 12 

f V 

t/2 




a p ) dt 


(A-13) 


and 


s - s Q + V • (^ k - a p ) t 


(A-14) 


However, recall that, in general. 


s = s 0 + f(t) 


(A-l 5) 


and, in the special case considered previously, 


s = s n + vt 


(A- 16) 


therefore, 


rt ,a a , 
v = V • (a k - a ) 
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The sine function for this case may be written as 


sinc [^T M ‘ k 2 P )] (A ' I7) 

Thus, the exposure at the hologram point, P, for this case becomes 

£(P) = m r | K c + E r E s sine K cos * J • (A-18) 

Looking again at the sine function of equation (A-17) and performing the indicated scalar 
products, it is seen that 


I - kr v 
sine 

A 

(K - s p) 

■ sinc [Ct 

- j (cos 7 - cos 6 ) 

/ J 

. (A- 19) 

L 2 


V /J 





in this form, the modifying function allows several facts to be noted by inspection. 
Except for the case of rv= 0 (i.e., stationary scene, because sine 0 = 1), the sine function 
has a maximum for only the following conditions: 

1 . That in which a^ is parallel to a p . 

2. That where V is perpendicular to the vector 

Consider that the vector a p always points toward the hologram plane (whatever its 
position); then, the condition that a^ and S p be parallel can be met only if ^ also points 

toward the hologram. It is seen that this condition is equivalent to the arrangement for 
the so-called direct- and diffuse-type holography. This maximum occurs at the point 
where the incident laser radiation intersects the hologram plane when projected through 
the scene point. It is noted that the sine function is a maximum for all orientations of 
the motion vector, V. 

The second condition for a maximum, that - a p ) always be perpendicular to 
the motion vector, V, is satisfied for all vectors (a^-ap) lying in a plane perpendicular 
to a y , since, for that condition, 
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£U> 


0 


(A-20) 


V • 



but, 


V • ^a^ - ap^ = v(cos 7 - cos 5) = 0 , (A-21) 

and this condition is satisfied whenever 5 = ± 7 . This condition is met by any vector, 

q, making an angle, 7 , with the vector, ^ r . The locus of these vectors, ap, forms a 

circular cone with a half angle equal to 7 about the motion vector and contains as an 
element (Fig. A-3). 


PLANE OF CONE BASEX TO V 



Figure A-3. Cone of constant fringe contrast. 
The zeros of the sine function, equation (A- 19), occur for 


sine 
n = 



± l, 2 ... 



sine n jr 


(A-22) 
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Therefore, the zeros occur for 


rkV 

__ 



(A-23) 


i.e., 


V • 



or 




nX 

VT 


(A- 24) 


(A-25) 


and 


nX 

COS 7 - cos 8 = 

117T v j. 


(A-26) 


or 


nX 

cos 5 n = COS 7 

V T 


(A-27) 


where 5 njr is that angle between vectors $ q and a y that produces a zero value for the 
sine function (i.e., that angle where no fringes will be recorded) 


From the above and Figure A-3, it is obvious that if a hologram is formed at any 
point in space, regions of the hologram having constant fringe contrast will form a conic 
section, since they represent the intersection of the hologram plane and conical loci 

mentioned above. If the vector Sp is allowed to take on various values, 6, in the range 
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7 K 6 < 7 + 2tt 


(A-28) 


for a constant with angle 7 , one can then construct all possible orientations 
of Sp about the motion vector, V. If for each specific value of 6 one allows it to rotate 
about the velocity vector, V, then one is able to construct all possible conical loci for ap. 
The results of this are shown in Figure A-4, where three possible values of 7 for aj, and 
all possible values of 5 for ^ for the specific 7 are taken. The values of the argument 
of the sine function for 7 = 0, n/4, and tt/ 2 are then plotted. 

Hologram-Oriented Coordinates . In the development of the previous subsection, 
the exposure was given by 


cos k [ s(t) - r] dt 


As before, the form of s(t) must be determined, but now as a function of the 
hologram. For this development, a fixed-coordinate system is used which has its origin at 
the center of the hologram and the z-axis normal to the plane of the hologram. This 
approach or development is the generalized derivation, since it may be used for the 
region near the hologram as well as the region far from the hologram by simply taking a 
sufficient number of terms in the series expansions that are developed. Its primary 
malady is that these very expansions tend to obscure the physical significance of the 
various steps. 

Note the diagram in Figure A-5 and let it be assumed that a hologram is to be 
formed in the x-y plane of this figure, which is a coordinate system with the origin at the 
center of the hologram. Analogous to the previous development, let 


(A-29) 


& (P) = m r 


K C + 


2E r E $ 


/ 

-t/2 


t . n (( 

s = s + S 


where s' is the distance from the scene point, P(x,y), of the hologram, and S" is the 
distance from the laser source to the scene point, S. (This is simply the total length of 
illumination called LQ + QP of the previous development.) Then, as before, because of 
the motion of S, the rate at which S" changes is 


d(S") 

dt 



V COS 7 
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Figure A-4. Variation of sine function argument with illumination direction. 



The coordinates for the point source which may be found from inspection of the 
diagram (Fig. A-5) are given by 

x g = r sin 6 cos , 

y s = r sin d sin ♦ , 

and (A-30) 

z„ = r cos d 
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Therefore, the distance, s', from the scene point, S, to a point, P(x,y), of the 
hologram is 


s' = (x - x s ) 2 + (y - y s ) 2 + z s 2 


( A-3 1 ) 


which becomes 


( y2 u2 9 y \ Vl 

— + — + 1 - — sin 6 cos 't' - — sin 8 sin 4') 
r 2 r 2 r r 


By expanding, one obtains 


»• = r(l - 


sin 6 cos - — sin 6 sin 4' 
r 


)♦£( 


1 - sin 2 0 cos 2 4' 


y 2 / 

+ — - (1 - sin 2 0 sin 2 4' 
2r 2 \ 


) -~ S1 


sin 2 8 sin 'I' cos 4' 


+ — - (sine cos 4' - sin 3 0 cos 3 ^ 
2r 2 V ) 


x 2 y /. 

— r- I si 

2r 2 ^ 


sin 8 sin 4* - sin 3 6 sin 4» cos 2 '!' 


xy 2 / 
2r 3 \ 


sin 0 cos 'l' - sin 3 8 sin 2 4' cos 4' 


+ — - | sin 8 sin 4' 

2r 2 V 


sin 3 # sin 3 '!' ) + (higher-order terms) . (A-32) 


If one assumes that x/r and y/r are sufficiently small so that only the first-order terms 
are needed, one obtains 


y . 

sin 8 cos 'I' - — sin 6 sm 
r 


(A-3 3) 
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To proceed, it becomes necessary to make some assumptions and then investigate 
some specific examples of the three orthogonal components of motion involved. 
Following are derivations of the exposures for linear motions of a constant phase laser 
source, $ , along each of three orthogonal directions. 

For a linear transverse motion (0), consider the point source, S, to be moving 
with a constant radius, tq, and in a constant 'I' plane, called 'J'q, such that 


e = 0 O + et 


then, 


sin 8 = sin (0g + 6 1) 


and 


sin 8 = sin 0 q cos 8 1 + cos 8 q sin 8 1 


However, 


(0t) 3 0 1) 5 

sin9t = »t + — 


and 


, (0t) 2 , (0t) 4 

cos 6 1 = 1 - — — + 

2 24 


If it is now assumed that 6 1 is small enough that only the first term in each series need 
be used, then one may write 


sin 8 = sin 8 q + Gt cos 0 q 
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then, 


sin 2 0 = sin 2 0Q + 2 dt sin 0q cos 0q 

and 

sin 3 0 = sin 3 0 q + 3 0 1 sin 2 0 q cos 0 q 

If these conditions are now substituted into equation (A-32) for the distance, s', 
one obtains 
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+ — (sin 0 o + 0 t cos 0 q) sin 'I' q 

2r 0 L 

- (sin 3 0Q + 3 0t sin 2 0g cos 0q) sin 'I 'q cos^'qJ 

+ |(sin 0 o + 0 1 cos 0 q) sin ^ 0 

- (sin 3 0 q + 3 0t sin 2 0 o cos 0q) sin 3 
+ (highcr-order terms) 

Again, only the first-order terms are taken, and equation (A-34) becomes 


X. 

= r 0 1 - — (sin 0Q + 0t cos 0q) cos 'I'q 

L r 0 


- — (sin 0 q + 0 $ t cos 0q) sin 'I'q 

r 0 


] 


or 

s’ = rQ - (x cos 'I'q + y sin 'I' q) sin 0 q - 0t (x cos 'I'q + y sin 'I' q) 

Now, for simplicity, define 

Sq' = tq - (x cos 'I'q + y sin 'I'q) sin 0q 

and one obtains 

s' = Sq' - 0 1 (x cos 'I'Q + y sin 'I'q) cos 00 


(A-34) 


cos 0 q . 


(A-35) 
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Then, the rate of change of phase between S and P, because of the motion in 
the 0 direction, is 





-kv d 


COS 0q 

(x cos •'I'o + y sin 'Fq) 


(A-36) 


where 


= 0r 

If one proceeds in a similar fashion, the derivations for small motions in 
the r and 'J' directions may be obtained. 

For linear radial motion, r, assume that the point source, S, moves with constant 
angles, 0 q and such that 


r = r n + rt 


then, 


in + rt 


or 


1 _ 1 
r r 0 


1 _ li + _ ^t) 3 

r o r o 2 r o 3 


and for rt < Tq , 


i-i<! 


2rt 3 Ft 2 

r o r o 2 
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If the motion, rt, is small compared with tq, one may again neglect higher-order 
terms, and by using the above expressions in equation (A-32), s r becomes 

s' = tq + rt - x sin 9 q cos 'k q - y sin 6 q sin 'k q 

[ x 2 y 2 

— (1 - sin 2 cos 2 'ko) + - — (1 - sin 2 sin 2 ^q) 

2r 0 2r 0 


xy 

— sin‘0r| sin 'tn cos 'kp 

r 0 


l r 0 


X 3 

+ — - (sin 9 q cos >k q - sin 3 (?Q cos 3 'kg) 

_ 2r 0 


x y , . 


+ - — - (sin 6 q sin ^0 " sin 3 0Q sin cos 2 ^q) 


XV 

+ - — ^ (sin 6 o cos 'I' g - sin 3 6 q sin 2 'k q cos 'kg) 


2r 0 2 


+ (sin 0 q sin 'kg - sin 3 0g sin 3 'kg 

2r 0~ 


1 - 


2rt 


Now, as before, taking only the first-order terms, one obtains 


s' = r 0 + rt - x sin 0 q cos 'kg - y sin 0 q sin ^q 


and, considering the direction of k, the phase rate is 


n r = - k 


l' r ' 

r 3s 

L dt 
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(A-37) 


n r = -kv r 


where 


V 


r 


For linear transverse motion, the third orthogonal component, it is assumed 
that the point source, S, moves with constant radius, tq, and angle, 0 q, such that 


— 'Fq + it 


then, 


sin ^ = sin 't'Q cos it + cos sin it 


and 


cos 'F = cos ^ q cos i t - sin 'F q sin i t 


Further, 


sin it 


it 


(it) 3 ^ (it) 5 

6 + "l20 


and 


cos it 


(it) 2 l (it) 4 

2 24 
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Again, for small motion, 'Pt < 'I'q 
sin ss 'i't 


and 


cos 'I' t « 1 


Using the proper substitution of these relations in equation (A-32), the following 
is obtained for s': 


s' = r« 1 - — sin 0g cos >Pg - — sin On sin 'Pg + — sin 0 q 4' t sin 'P q 
1 tq rn rn 


r 0 


y • x* 

- — sin 8 q 'P t cos ^ ft + r — — f 1 - sin 2 8 n cos 2 'P 


2r, 


7(“ 


' 0 COS T 0 


■) 


+ — — - sin 2 0 a sin^Q cos 'I 'q + 

2tq 2r, 


H 1 - 


sin 2 0 q sin 2 'P 


’°) 


y 2 xy 

— sin 2 0g 'P t sin 'P q cos 'P q - — sin 2 0 q simPg cos 'Pg 
r r Q 


- —■ sin 2 0g 'Pt (cos 2 'Pg - sin 2 'Pg ) 

r 0 V / 


o V 


+ - — - (sin 0 q cos 'Pq - 'Pt sin 0 q sin <Pq - sin 3 0gcos 'Pg 
2r r 


+ 3 'Pt sin 3 'Pg cos 2 'P 


°) + ^[ Sin0 ° Sin ' I, ° 


'Pt sin 0 q cos 'Pq 


- sin 3 0g sin 'Pg cos 2 'Pg + 'Pt sin 3 0g (2 sin 2 'Pg cos 'Pg - cos 3 'P 


•] 
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xy 

+ sin 0 q cos 'I' q - 'kt sin 0g sin +g - sin 3 0g sin 2 ^g cos 'kg 

2 r o 3 L 


- sin 3 0Q (2 sin cos 2 


*0 “ sin 3 'kg) J 


+ (sin 0 q sin ^ q + >kt sin 0q cos 'kg - sin 3 0g sin 3 4g 

2rn 


- 3 'kt sin 3 0q sin 2 'I'g cos 'I'g) 


Again, taking only first-order terms, one obtains 


s' = rg - x sin 0 g cos 'I' g - y sin 0 sin 'I' g + sin 0 g 'i' t sin 'I' g 


y sin 0 g ^ t cos 'I' g 


and again noting the direction of k, 


£2 T = - k V T 
^ 'k 


(x sin ^g - y cos 'kg) 
r 0 


(A-38) 


where 


v = 'k rg sin 0 g 


Any general linear motion may be resolved into the three orthogonal components 
used in the above development; i.e., (V p V^). The constant portion will be a sum of 

the separate contributions to the phase shifts, and the time-variant portions will be 
formed by summing the time variant phases found for the three orthogonal components. 
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The recorded fringe amplitude will then be modified by a sine function whose argument 


— - (£2 r + n e + £2^)- 


Therefore, combining the time-variant and the constant portions of the phase 
shifts, the total phase retardation from the laser to the hologram by way of the scene is 


ks(t) - k |^Sq" + ^ V - a^r J + k sq' + (£2 r + £2 g + £2^)' 
ks(t) = kjjSg' + Sq"J + (kv cos 7 + £2)r , 


ks(t) = k sg + (kv cos 7 + £2 )t 


where 


£2 = £ 2 r + £ 2 ^ + £ 2 ^ 


(A- 39 ) 


The exposure for this case of linear motion may be found by analogy with the 
exposure for the linear case represented by equation (A- 7 ), where, by analogy, it can be 
seen that 


kv -*■ kv cos 7 + £2 


(A- 40 ) 


The exposure for this hologram-oriented coordinate case then becomes 


& (P) = mr + E r E s sine — (v cos y + £2 ) 


)j COS <^| 


(A- 41 ) 


where the constants, m and Kg,, are as defined earlier. Then the fringes are modulated by 
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sine 


k r 

T 


(v cos 7 + £2) 


■] 


kr 


sine 


[■ 


v cos y 


V 0 

— (x cos 'I'q + y sin 'i'Q) cos 0 q 


- V 


r 


(x sin 'I'q - y cos 'I'q) 


(A-42) 


This is then the form of the cosine fringe modulation function in terms of the 
hologram-oriented coordinates. 


If one lets r -*■ the sine function of equation (A-42) becomes 


sine 


[ 


kr 

— (v cos 7 
2 


+ S2 ) 

( 

j 


sine 


(v cos 7 - V 


- 

r> 


(A-43) 


Recall that 8 is the angle between the velocity vector and a unit vector in the 
radial direction of the previous vector analysis development; then, by analogy, 


V r = V ■ ap = v cos 5 


(A-44) 


and equation (A-43) becomes 


sine 


kr 


(v cos y + £2 ) 


I 


krv 

sine 

2 


(cos y - COS 8) 


(A45) 


and one has the expected result. 


Resultant Effects of Linear Motion on the Reconstructed Wavefront 

The previous section complements this section in showing that the effect of 
motion of a point source scene during the formation of a hologram is a modification of 
the recorded fringe contrast at various points of the hologram. Furthermore, the effects 
on the hologram’s recorded fringes caused by the scene motion may be derived by 
considering the motion of only a single point of a rigid moving scene. 
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Now the characteristics of the image reconstructed from a hologram having such 
motion-modified fringes will be investigated. A comparison of this image with the original 
object will give a point-spread function for the motion holography case. The approach 
followed will be that of the currently existing method presented in Reference 12. 

Consider Figure A-6, where a particular source distribution in the vicinity of the 
scene is assumed. It is desirable then to find the resulting amplitude and phase 

distribution on the hologram at point Q. If it can be shown that this is identical to the 

amplitude and phase distribution of the reconstructed field, it will prove that the 
assumed source will describe the reconstructed virtual image. One chooses a uniform line 
source lying along the motion vector as the assumed source (Fig. A-6), i.e., the Z-axis in 
this diagram. The uniform line source of the length, AC, is centered at the t = 0 scene 
point, sq, and has a linear phase variation, kbz, where b is some proportionality constant. 

The resulting field will be symmetrical with respect to revolution about the Z-axis. Then, 
it is desired to know the field at a point, Q, on the hologram^ a distance, qp » AJ2, 

from the source and at an angle, 5, from the Z-axis. 



Figure A-6. Geometry for reconstruction analysis. 
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We consider that the field dE at the point, Q, is due to the contributions of point 

sources in the infinitesimal length, dZ, at a distance, Z, from the origin. Then, dE is given 
by 


A 

dE = — exp [i(co t + kbZ - kq )] dZ 
^0 

From the diagram, 

qg = q + Z cos 6 ; 


(A-46) 


(A-47) 


thus. 


q = qg - Z cos 6 


Then, 


dE - — exp [i(cj t + kbZ - kq n + kZ cos 5)] dZ , (A-48) 

qo 

where A is a constant involving amplitude. 

Now, by integrating, one finds the total field at Q 


E(Q) = — exp i (cot - 

% 


AC/2 

kq Q ) f exp [ik (b 
AC/2 


- cos 6) Z] dZ . (A-49) 


Applying Euler’s formula to the integral and setting in the limits of integration, 


AS 

^ . 2 sin [k (b - cos 6) — 

E(Q) = — exn i(co t - kq n )'l ^ 

U I k (b - cos 6) 


(A-50) 
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Upon multiplying the brackets by AC/AC = 1, 


E(Q) = — exp i (cut - kqQ) AC sine k(b - cos 8) — 

Qo L 2 


AC' 

y 


(A-5 1 ) 


Therefore, E(Q) represents a homocentric wave, centered at the origin point, 0, with an 
amplitude weighting given by 


sine 


kAC 


(b - cos 8) 


(A-52) 


Equation (A-5 1 ) is the description of the field at the point, Q, of the hologram 
plane resulting from the assumed source of a uniform line. 

It is now appropriate to determine the amplitude and phase distribution of the 
hologram’s reconstructed field and then compare the reconstructed (virtual image) field 
with that of the assumed source. If the two are identical, it will be obvious that the 
assumed source will describe the reconstructed virtual image. 

For the reconstructed hologram, the amplitude transmission factor may be written 
as 


T a = T 0 - M 


(A-53) 


The exposure, £ , is given by equation (A-l 8), 


£ = m t I 

iK^ + E r E s sine 

krV 

y 

(\ - %) 

cos (0 r - 0 S ) 


[ 


- 



or 


£ 


Kl 


+ k 2 


krV 





cos (0 r - <P S ) 


(A-54) 
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Then, from the amplitude transmission factor, one can write 



Now, since equations (A-53), (A-54), and (A-55) describe the exposure and 
development of the field resulting from a motion scene, the reconstructed virtual image 
field (first-order diffraction) has an amplitude proportional to the magnitude of the 
spatial variations in the amplitude transmission factor and is, therefore, proportional to 


krV 




kvr 

sine — (cos 7 - cos 8 ) 


(A-56) 


Therefore, since the fringe phase (<I> r *<t> $ ) is independent of the motion, V, the 

hologram will reconstruct a homocentric wave with the center at the t = 0 position of the 
scene point. 

Then, comparing equations (A-52) and (A-56), one sees that, if A£ = vt and 
b = cos 7 , the phase fronts and amplitude distributions are identical for the assumed 
source and the reconstructed (virtual image) field of the hologram. Therefore, the virtual 
image will be a line source of length, vt, and with a phase equal to KZ cos 7, where Z is 
the distance from the center of the source. In other words, the image will have the same 
blur as results from any detector (such as a conventional photograph) having an 
integrating time of r. Furthermore, the phase relationships are preserved between the 
impulses that sequentially formed the hologram. 
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